Investigating distributions of inhaled aerosols in the lungs of post-COVID-19 clusters through a unified imaging and modeling approach.

BACKGROUND: Recent studies, based on clinical data, have identified sex and age as significant factors associated with an increased risk of long COVID. These two factors align with the two post-COVID-19 clusters identified by a deep learning algorithm in computed tomography (CT) lung scans: Cluster 1 (C1), comprising predominantly females with small airway diseases, and Cluster 2 (C2), characterized by older individuals with fibrotic-like patterns. This study aims to assess the distributions of inhaled aerosols in these clusters. METHODS: 140 COVID survivors examined around 112 days post-diagnosis, along with 105 uninfected, non-smoking healthy controls, were studied. Their demographic data and CT scans at full inspiration and expiration were analyzed using a combined imaging and modeling approach. A subject-specific CT-based computational model analysis was utilized to predict airway resistance and particle deposition among C1 and C2 subjects. The cluster-specific structure and function relationships were explored. RESULTS: In C1 subjects, distinctive features included airway narrowing, a reduced homothety ratio of daughter over parent branch diameter, and increased airway resistance. Airway resistance was concentrated in the distal region, with a higher fraction of particle deposition in the proximal airways. On the other hand, C2 subjects exhibited airway dilation, an increased homothety ratio, reduced airway resistance, and a shift of resistance concentration towards the proximal region, allowing for deeper particle penetration into the lungs. CONCLUSIONS: This study revealed unique mechanistic phenotypes of airway resistance and particle deposition in the two post-COVID-19 clusters. The implications of these findings for inhaled drug delivery effectiveness and susceptibility to air pollutants were explored.

Transport and deposition of beclomethasone dipropionate drug aerosols with varying ethanol concentration in severe asthmatic subjects.

This study aims to assess the effects of varying an ethanol co-solvent on the deposition of drug particles in severe asthmatic subjects with distinct airway structures and lung functions using computational fluid dynamics. The subjects were selected from two quantitative computed tomography imaging-based severe asthmatic clusters, differentiated by airway constriction in the left lower lobe. Drug aerosols were assumed to be generated from a pressurized metered-dose inhaler (MDI). The aerosolized droplet sizes were varied by increasing the ethanol co-solvent concentration in the MDI solution. The MDI formulation consists of 1,1,2,2-tetrafluoroethane (HFA-134a), ethanol, and beclomethasone dipropionate (BDP) as the active pharmaceutical ingredient. Since HFA-134a and ethanol are volatile, both substances evaporate rapidly under ambient conditions and trigger condensation of water vapor, increasing the size of aerosols that are predominantly composed of water and BDP. The average deposition fraction in intra-thoracic airways for severe asthmatic subjects with (or without) airway constriction increased from 37%±12 to 53.2%±9.4 (or from 20.7%± 4.6 to 34.7%±6.6) when the ethanol concentration was increased from 1 to 10%wt/wt. However, when the ethanol concentration was further increased from 10 to 20%wt/wt, the deposition fraction decreased. This indicates the importance of selecting appropriate co-solvent amounts during drug formulation development for the treatment of patients with narrowed airway disease. For severe asthmatic subjects with airway narrowing, the inhaled aerosol may benefit from a low hygroscopic effect by reducing ethanol concentration to penetrate the peripheral region effectively. These results could potentially inform the selection of co-solvent amounts for inhalation therapies in a cluster-specific manner.

Characterizing Subjects Exposed to Humidifier Disinfectants Using Computed-Tomography-Based Latent Traits: A Deep Learning Approach.

Around nine million people have been exposed to toxic humidifier disinfectants (HDs) in Korea. HD exposure may lead to HD-associated lung injuries (HDLI). However, many people who have claimed that they experienced HD exposure were not diagnosed with HDLI but still felt discomfort, possibly due to the unknown effects of HD. Therefore, this study examined HD-exposed subjects with normal-appearing lungs, as well as unexposed subjects, in clusters (subgroups) with distinct characteristics, classified by deep-learning-derived computed-tomography (CT)-based tissue pattern latent traits. Among the major clusters, cluster 0 (C0) and cluster 5 (C5) were dominated by HD-exposed and unexposed subjects, respectively. C0 was characterized by features attributable to lung inflammation or fibrosis in contrast with C5. The computational fluid and particle dynamics (CFPD) analysis suggested that the smaller airway sizes observed in the C0 subjects led to greater airway resistance and particle deposition in the airways. Accordingly, women appeared more vulnerable to HD-associated lung abnormalities than men.

A computed tomography imaging-based subject-specific whole-lung deposition model.

The respiratory tract is an important route for beneficial drug aerosol or harmful particulate matter to enter the body. To assess the therapeutic response or disease risk, whole-lung deposition models have been developed, but were limited by compartment, symmetry or stochastic approaches. In this work, we proposed an imaging-based subject-specific whole-lung deposition model. The geometries of airways and lobes were segmented from computed tomography (CT) lung images at total lung capacity (TLC), and the regional air-volume changes were calculated by registering CT images at TLC and functional residual capacity (FRC). The geometries were used to create the structure of entire subject-specific conducting airways and acinar units. The air-volume changes were used to estimate the function of subject-specific ventilation distributions among acinar units and regulate flow rates in respiratory airway models. With the airway dimensions rescaled to a desired lung volume and the airflow field simulated by a computational fluid dynamics model, particle deposition fractions were calculated using deposition probability formulae adjusted with an enhancement factor to account for the effects of secondary flow and airway geometry in proximal airways. The proposed model was validated in silico against existing whole-lung deposition models, three-dimensional (3D) computational fluid and particle dynamics (CFPD) for an acinar unit, and 3D CFPD deep lung model comprising conducting and respiratory regions. The model was further validated in vivo against the lobar particle distribution and the coefficient of variation of particle distribution obtained from CT and single-photon emission computed tomography (SPECT) images, showing good agreement. Subject-specific airway structure increased the deposition fraction of 10.0-µm particles and 0.01-µm particles by approximately 10%. An enhancement factor increased the overall deposition fractions, especially for particle sizes between 0.1 and 1.0 µm.

Longitudinal Imaging-Based Clusters in Former Smokers of the COPD Cohort Associate with Clinical Characteristics: The SubPopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS).

PURPOSE: Quantitative computed tomography (qCT) imaging-based cluster analysis identified clinically meaningful COPD former-smoker subgroups (clusters) based on cross-sectional data. We aimed to identify progression clusters for former smokers using longitudinal data. PATIENTS AND METHODS: We selected 472 former smokers from SPIROMICS with a baseline visit and a one-year follow-up visit. A total of 150 qCT imaging-based variables, comprising 75 variables at baseline and their corresponding progression rates, were derived from the respective inspiration and expiration scans of the two visits. The COPD progression clusters identified were then associated with subject demography, clinical variables and biomarkers. RESULTS: COPD severities at baseline increased with increasing cluster number. Cluster 1 patients were an obese subgroup with rapid progression of functional small airway disease percentage (fSAD%) and emphysema percentage (Emph%). Cluster 2 exhibited a decrease of fSAD% and Emph%, an increase of tissue fraction at total lung capacity and airway narrowing over one year. Cluster 3 showed rapid expansion of Emph% and an attenuation of fSAD%. Cluster 4 demonstrated severe emphysema and fSAD and significant structural alterations at baseline with rapid progression of fSAD% over one year. Subjects with different progression patterns in the same cross-sectional cluster were identified by longitudinal clustering. CONCLUSION: qCT imaging-based metrics at two visits for former smokers allow for the derivation of four statistically stable clusters associated with unique progression patterns and clinical characteristics. Use of baseline variables and their progression rates enables identification of longitudinal clusters, resulting in a refinement of cross-sectional clusters.

Assessment and Validation of a Hygroscopic Growth Model with Different Water Activity Estimation Methods.

Hygroscopic growth models are currently of interest as aids for targeting the deposition of inhaled drug particles in preferred areas of the lung that will maximize their pharmaceutical effect. Mathematical models derived to estimate hygroscopic growth over time have been previously developed but have not been thoroughly validated. For this study, model validation involved a comparison of modeled values to measured values when the growing droplet had reached equilibrium. A second validation process utilized a novel system to measure the growth of a droplet on a microscope coverslip relative to modeled values when the droplet is undergoing the initial rapid growth phase. Various methods currently used to estimate the water activity of the growing droplet, which influences the droplet growth rate, were also compared. Results indicated that a form of the hygroscopic growth model that utilizes coupled-differential equations to estimate droplet diameter and temperature over time was valid throughout droplet growth until it reached its equilibrium size. Accuracy was enhanced with the use of a polynomial expression to estimate water activity relative to the use of a simplified estimate of water activity based on Raoult's Law. Model accuracy was also improved when constraining the film of salt solution surrounding the dissolving salt core at saturation.

An integrated 1D breathing lung simulation with relative hysteresis of airway structure and regional pressure for healthy and asthmatic human lungs.

This study aims to develop a one-dimensional (1D) computational fluid dynamics (CFD) model with dynamic airway geometry that considers airway wall compliance and acinar dynamics. The proposed 1D model evaluates the pressure distribution and the hysteresis between the pressure and tidal volume (Vtidal) in the central and terminal airways for healthy and asthmatic subjects. Four-dimensional CT images were captured at 11-14 time points during the breathing cycle. The airway diameter and length were reconstructed using a volume-filling method and a stochastic model at respective time points. The obtained values for the airway diameter and length were interpolated via the Akima spline to avoid unboundedness. A 1D energy balance equation considering the effects of wall compliance and parenchymal inertance was solved using the efficient aggregation-based algebraic multigrid solver, a sparse matrix solver, reducing the computational costs by around 90% when compared with the generalized minimal residual solver. In the Vtidal versus displacement in the basal direction (z-coordinate), the inspiration curve was lower than the expiration curve, leading to relative hysteresis. The dynamic deformation model was the major factor influencing the difference in the workload in the central and terminal airways. In contrast, wall compliance and parenchymal inertance appeared only marginally to affect the pressure and workload. The integrated 1D model mimicked dynamic deformation by predicting airway diameter and length at each time point, describing the effects of wall compliance and parenchymal inertance. This computationally efficient model could be utilized to assess breathing mechanism as an alternative to pulmonary function tests.NEW & NOTEWORTHY This study introduces a one-dimensional (1D) computational fluid dynamics (CFD) model mimicking the realistic changes in diameter and length in whole airways and reveals differences in lung deformation between healthy and asthmatic subjects. Utilizing computational models, the effects of parenchymal inertance and airway wall compliance are investigated by changing ventilation frequency and airway wall elastance, respectively.

Transport and deposition of hygroscopic particles in asthmatic subjects with and without airway narrowing.

This study numerically investigates the effect of hygroscopicity on transport and deposition of particles in severe asthmatic lungs with distinct airway structures. The study human subjects were selected from two imaging-based severe asthmatic clusters with one characterized by non-constricted airways and the other by constricted airways in the lower left lobe (LLL). We compared the deposition fractions of sodium chloride (NaCl) particles with a range of aerodynamic diameters (1-8 µm) in cluster archetypes under conditions with and without hygroscopic growth. The temperature and water vapor distributions in the airways were simulated with an airway wall boundary condition that accounts for variable temperature and water vapor evaporation at the interface between the lumen and the airway surface liquid layer. On average, the deposition fraction increased by about 6% due to hygroscopic particle growth in the cluster subjects with constricted airways, while it increased by only about 0.5% in those with non-constricted airways. The effect of particle growth was most significant for particles with an initial diameter of 2 µm in the cluster subjects with constricted airways. The effect diminished with increasing particle size, especially for particles with an initial diameter larger than 4 µm. This suggests the necessity to differentiate asthmatic subjects by cluster in engineering the aerosol size for tailored treatment. Specifically, the treatment of severe asthmatic subjects who have constricted airways with inhalation aerosols may need submicron-sized hygroscopic particles to compensate for particle growth, if one targets for delivering to the peripheral region. These results could potentially inform the choice of particle size for inhalational drug delivery in a cluster-specific manner.

Combining existing numerical models with data assimilation using weighted least-squares finite element methods.

A new approach has been developed for combining and enhancing the results from an existing computational fluid dynamics model with experimental data using the weighted least-squares finite element method (WLSFEM). Development of the approach was motivated by the existence of both limited experimental blood velocity in the left ventricle and inexact numerical models of the same flow. Limitations of the experimental data include measurement noise and having data only along a two-dimensional plane. Most numerical modeling approaches do not provide the flexibility to assimilate noisy experimental data. We previously developed an approach that could assimilate experimental data into the process of numerically solving the Navier-Stokes equations, but the approach was limited because it required the use of specific finite element methods for solving all model equations and did not support alternative numerical approximation methods. The new approach presented here allows virtually any numerical method to be used for approximately solving the Navier-Stokes equations, and then the WLSFEM is used to combine the experimental data with the numerical solution of the model equations in a final step. The approach dynamically adjusts the influence of the experimental data on the numerical solution so that more accurate data are more closely matched by the final solution and less accurate data are not closely matched. The new approach is demonstrated on different test problems and provides significantly reduced computational costs compared with many previous methods for data assimilation. Copyright © 2016 John Wiley & Sons, Ltd.

Left ventricular flow analysis: recent advances in numerical methods and applications in cardiac ultrasound.

The left ventricle (LV) pumps oxygenated blood from the lungs to the rest of the body through systemic circulation. The efficiency of such a pumping function is dependent on blood flow within the LV chamber. It is therefore crucial to accurately characterize LV hemodynamics. Improved understanding of LV hemodynamics is expected to provide important clinical diagnostic and prognostic information. We review the recent advances in numerical and experimental methods for characterizing LV flows and focus on analysis of intraventricular flow fields by echocardiographic particle image velocimetry (echo-PIV), due to its potential for broad and practical utility. Future research directions to advance patient-specific LV simulations include development of methods capable of resolving heart valves, higher temporal resolution, automated generation of three-dimensional (3D) geometry, and incorporating actual flow measurements into the numerical solution of the 3D cardiovascular fluid dynamics.

A model of filiform hair distribution on the cricket cercus.

Crickets and other orthopteran insects sense air currents with a pair of abdominal appendages resembling antennae, called cerci. Each cercus in the common house cricket Acheta domesticus is covered with between 500 to 750 filiform mechanosensory hairs. The distribution of the hairs on the cerci, as well as the global patterns of their movement axes, are very stereotypical across different animals in this species, and the development of this system has been studied extensively. Although hypotheses regarding the mechanisms underlying pattern development of the hair array have been proposed in previous studies, no quantitative modeling studies have been published that test these hypotheses. We demonstrate that several aspects of the global pattern of mechanosensory hairs can be predicted with considerable accuracy using a simple model based on two independent morphogen systems. One system constrains inter-hair spacing, and the second system determines the directional movement axes of the hairs.